Mapa de plagas

Column

Filters

Column

Datatable

Interactive map

Mapa especifico

Row

mapa

---
title: "Mapa interactivo"
author: "Gorky Florez Castillo"
date: "31/7/2021"
output:  
  flexdashboard::flex_dashboard:
    theme: paper
    source_code: embed
---
```{r setup, include=FALSE}
# setting up R Markdown options

# We want to hide the code and only see the results
knitr::opts_chunk$set(echo = F)

# We don't want to see any warnings from our code
knitr::opts_chunk$set(warning = F)

# We don't want to see any messages
knitr::opts_chunk$set(message = F)

```

```{r install_packages}
# You must have the flexdashboard package installed
# Before knitting this R Markdown file
# install.packages("flexdashboard")

# This function checks if you don't have the correct packages installed yet
# If not, it will install it for you
packages <- c("tidyverse", "flexdashboard",
              "crosstalk", "leaflet", "DT")

if (length(setdiff(packages, rownames(installed.packages()))) > 0) {
  install.packages(setdiff(packages, rownames(installed.packages())), repos = "http://cran.us.r-project.org")  
}
library(tidyverse)
library(flexdashboard)
library(crosstalk)  
library(leaflet)   
library(DT)
library(leaflet)
library(dplyr)
library(leaflet.extras)
library(htmltools)
library(htmlwidgets)
library(crosstalk)
library(rgdal)
library(openxlsx)
library(png)  
library(broom)
library(tidyverse)
library(leaflet.providers)
library(leafem)
library(sp)
library(sf)#Manejo de informacion espacial
library(mapview)#Para visualizacion de datos espaciales
library(RColorBrewer) #Paleta de Colores
library(viridis)
library(Rgb)
library(ggplot2)#Para distintos graficos incluso mapas
library(raster)#Para leer archivos raster
library(rmarkdown)

```

```{r cargar_limpiar_datos}
# Este ejemplo se centra en Luisiana
# Vuelve a esto más tarde si quieres
# para cambiar el estado que estamos visualizando
# NOTA: Deberá cambiar el inicio
# latitud y longitud del código del mapa a continuación
data <- read.xlsx("Data excel/full_data.xlsx", sheet="Sheet1")  %>% 
  filter(Lote=="LA") 

st <- SharedData$new(data)
```

Mapa de plagas {data-icon="ion-stats-bars"}
=====================================  

Column {data-width=150}
-------------------------------------

### Filters

```{r filter_section}
filter_select(
  id = "Lote_Sect",
  label = "Lotes",
  sharedData = st,
  group = ~Lote_Sect
)

bscols(
  filter_checkbox(
    id = "Fecha_muestreo",
    label = "Fechas",
    sharedData = st,
    group = ~Fecha_muestreo))

bscols(
  filter_slider(
    id = "Porcentaje_Botrytis",
    label = "Porcentaje de la plaga",
    sharedData = st,
    column = ~Porcentaje_Botrytis,
    step = 10,
    round = TRUE,
    sep = "",
    ticks = FALSE))
```


Column {data-width=500}
-------------------------------------

### Datatable

```{r filterable_table}
datatable(st)
```

### Interactive map

```{r interactive_map}
#Definiendo el Logo
m="https://images.vexels.com/media/users/3/143561/isolated/preview/afa3aa927b63061e3b0222b7dab9cdbf-ubicaci--n-n--utica-norte-flecha-vintage-by-vexels.png"


st %>% 
  leaflet() %>%
   addLogo(m,url = "https://images.vexels.com/media/users/3/143561/isolated/preview/afa3aa927b63061e3b0222b7dab9cdbf-ubicaci--n-n--utica-norte-flecha-vintage-by-vexels.png",
          position = "topleft",
          offset.x = 50,
          offset.y = 10,
          width = 100,
          height = 100)%>%
  addControl(html = "

Mapa ataque de plagas en cultivos

", position = "topright")%>% addDrawToolbar(targetGroup = "Graficos", editOptions = editToolbarOptions(selectedPathOptions = selectedPathOptions()))%>% addMeasure(position = "topleft", primaryLengthUnit = "meters", primaryAreaUnit = "sqmeters", activeColor = "#3D535D", completedColor = "#7D4479")%>% addLayersControl(baseGroups = c("Satellite", "OSM", "OTM", "ENGWM", "CDBP"), position = "topright", options = layersControlOptions(collapsed = T))%>% addProviderTiles(providers$Esri.WorldImagery, group = "Satellite")%>% addProviderTiles(providers$OpenStreetMap, group = "OSM")%>% addProviderTiles(providers$OpenTopoMap, group = "OTM")%>% addProviderTiles(providers$CartoDB.Positron, group = "CDBP") %>% addProviderTiles(providers$Esri.NatGeoWorldMap, group = "ENGWM")%>% setView(-77.6850487, -10.716505, zoom = 13) %>% addCircles(popup = ~paste0("", data$Lote_Sect ,"
", data$Fecha_muestreo, "
", data$Porcentaje_Botrytis), radius = data$Porcentaje_Botrytis, stroke=FALSE, opacity=.5)%>% addScaleBar(position = "bottomright",options = scaleBarOptions(maxWidth = 100, metric = TRUE, imperial = TRUE, updateWhenIdle = TRUE))%>% addMiniMap(tiles = providers$Esri.WorldImagery, toggleDisplay = TRUE) ``` Mapa especifico {data-icon="ion-stats-bars"} ===================================== Row {data-width=150} ------------------------------------- ### mapa ```{r interacti_map} #Definiendo el Logo m="https://images.vexels.com/media/users/3/143561/isolated/preview/afa3aa927b63061e3b0222b7dab9cdbf-ubicaci--n-n--utica-norte-flecha-vintage-by-vexels.png" st %>% leaflet() %>% addLogo(m,url = "https://images.vexels.com/media/users/3/143561/isolated/preview/afa3aa927b63061e3b0222b7dab9cdbf-ubicaci--n-n--utica-norte-flecha-vintage-by-vexels.png", position = "topleft", offset.x = 50, offset.y = 10, width = 100, height = 100)%>% addControl(html = "

Mapa ataque de plagas en cultivos

", position = "topright")%>% addDrawToolbar(targetGroup = "Graficos", editOptions = editToolbarOptions(selectedPathOptions = selectedPathOptions()))%>% addMeasure(position = "topleft", primaryLengthUnit = "meters", primaryAreaUnit = "sqmeters", activeColor = "#3D535D", completedColor = "#7D4479")%>% addLayersControl(baseGroups = c("Satellite", "OSM", "OTM", "ENGWM", "CDBP"), position = "topright", options = layersControlOptions(collapsed = T))%>% addProviderTiles(providers$Esri.WorldImagery, group = "Satellite")%>% addProviderTiles(providers$OpenStreetMap, group = "OSM")%>% addProviderTiles(providers$OpenTopoMap, group = "OTM")%>% addProviderTiles(providers$CartoDB.Positron, group = "CDBP") %>% addProviderTiles(providers$Esri.NatGeoWorldMap, group = "ENGWM")%>% setView(-77.6850487, -10.716505, zoom = 13) %>% addCircles(popup = ~paste0("", data$Lote_Sect ,"
", data$Fecha_muestreo, "
", data$Porcentaje_Botrytis), radius = data$Porcentaje_Botrytis, stroke=FALSE, opacity=.5)%>% addScaleBar(position = "bottomright",options = scaleBarOptions(maxWidth = 100, metric = TRUE, imperial = TRUE, updateWhenIdle = TRUE))%>% addMiniMap(tiles = providers$Esri.WorldImagery, toggleDisplay = TRUE) ```